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MODEL DEVELOPMENT OF SUPERSONIC TROUGH 
WIND WITH SHOCKS 

J. M. Grebowsky 
ABSTRACT 

The time dependent one dimensional hydrodynamic equations describe the 
evolution of the thermal plasma flow along closed magnetic field lines outside 
of the plasmasphere. The convection of the supersonic polar wind onto a closed 
field line results in the assumed formation of collisionless plasma shocks. These 
shocks move earthward as the field line with its "frozen-in’' plasma remains 
fixed or contracts with time to smaller L coordinates. The high equatorial 
plasma temperature (of the order of electron volts) produced by the shock process 
decreases with time if the flow is isothermal but it will increase if the contraction 
is under adiabatic conditions. Assuming adiabaticity a peak in the temperature 
forms at the equator in conjunction with a depression in the ion density. After 
an initial contraction, if the flux tube drifts to higher L coordinates the direction 
of the shock motion can be reversed so that the supersonic region wiH expand 
along the field line towards the state characterizing the supersonic polar wind. 
A rapid expansion will lower the equatorial density while the temperature de- 
creases with time under adiabatic but not isothermal conditions. 
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MODEL DEVELOPMENT OF SUPERSONIC TROUGH 


WIND WITH SHOCKS 


Introduction 

It is well established that the plasmasphere abruptly terminates at a mag- 
netic field aligned shell which intersects the equatorial plane at a geocentric 
distance of a few earth radii. This boundary, the plasmapause, corresponding 
to a rapid decrease of the light ion density with increasing L coordinate also 
appears to separate two regions of distinctly different light ion flow behavior 
(Mayr et al., 1970; Brinton et al., 1971, Banks et al., 1971). Within the plasma- 
sphere the macroscopic flow of thermal H + ions along the field lines is subsonic 
with the direction of the flow either away from or toward the earth depending on 
the phase of the diurnal rotation. On the other hand, the closed field lines out- 
side of the plasmasphere (i.e., the field lines associated with the light ion trough) 
are characterized by regions of upward directed flows of the thermal H + ions 
which approach or even exceed the sound speed. 

The upward directed supersonic H + fluxes along the trough field lines are analo- 
gous to and are directly related to the supersonic polar wind flowing along open polar 
lines (Banks and Holzer, 1969a). The relation between the trough wind and the polar 
wind is depicted in Figure 1. Through field line reconnection the supersonically flow- 
ing polar wind convects onto closed magnetic field lines extending to the earth's 
magnetospheric tail. The upward directed supersonic H + fluxes from the northern 
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and southern hemispheres will initially merge near the magnetic equator. In- 
stabilities will prevent the interpretation of these oppositely directed streams 
(Koons and Schulz, 1970; Schulz and Koons, 1972) converting the ordered flow 
motion into thermal energy. 

If the transition from supersonic to subsonic flow along a trough field line 
occurs rapidly enough to be treated as a shock discontinuity (Banks et al., 1971) 
two collisionless shocks - one in each hemisphere - will be formed initially on 
the closed field line resulting from reconnection in the magnetosphere tail. 

These shocks will propagate away from the equator along the field line towards 
the earth as the ambient plasma, frozen to the magnetic field lines, drifts 
towards the dayside magnetopause where the closed field line opens via connec- 
tion to the interplanetary field resulting in the expansive polar wind. 

Banks et al. (1971) have computed the motion of these shocks along a mag- 
netic dipole field line which is fixed in space — i.e., along a field line character- 
ized by a constant L coordinate. These computations yielded the important 
result that for an isothermal supersonic H + flow along a fixed L field line the 
characteristic travel time of the shocks from the equator to the base of the 
protonosphere can be much larger than the average time required for a closed 
flux tube with its associated plasma to convect from the magnetosphere tail to the 
dayside magnetopause. Hence within the light ion trough the lower protosphere 
should on the average be characterized by an upward directed supersonic flow of 
the H + ions. 
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This paper will expand upon the study of Banks et al. (1971) by explicitly 
evaluating the protonospheric flow development along a closed trough field line for 
both isothermal and adiabatic flow states. Also since the ambient plasma can 
drift perpendicular to the direction of the magnetic field, the compression or 
expansion of a magnetic flux tube and its associated plasma as occurs during a 
convection cycle must be considered in a complete study of the trough flow. The 
effects of the cross L drift on the supersonic flow of H + and on the shock motion 
will be explored. 

Theory 

The hydrodynamic description has been shown by Holzer et al. (1971) to be 
a good approximation to the supersonic flow state characterizing the polar wind. 
Such a description was used in previous studies of the polar wind (Banks and 
Holzer, 1969a, 1969b; Banks et al., 1971) and will be used to approximate the 
supersonic flow region on trough field lines. The subsonic flow downstream of 
the shocks has been randomized by the shock process and will be characterized 
by an isotropic pressure so that the flow everywhere along the trough field line 
is describable by the hydrodynamic relations. 

Previous model studies of the supersonic flow of the polar and trough winds 
used the steady state one dimensional form of the hydrodynamic equations to 
describe the flow along a field line. This is valid when the L coordinate is 
constant and when steady state conditions prevail at the base of the field line 
because any departure of the supersonic region from static equilibrium results 
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in a rapid relaxation to the steady state configuration. However, fixing the frame 
of reference in the plane containing the trough field line under investigation, the 
field aligned subsonic flow associated with the shock transition requires the use 
of the time dependent hydrodynamic equations since the shock position varies 
with time. Hence attention will be devoted to solving this form of the flow equa- 
tions. The one dimensional time dependent equations can also be extended to 
study the flow variations produced by cross L drifts — this is to be considered 
in a later section. 

Along a field line (assumed dipolar) characterized by an L coordinate fixed 
in time, the hydrodynamic equations can be resolved into components along the 
field direction with two independent variables; the time t and the distance s 
along the field line from an arbitrary base level. Considering the protonosphere 
where chemical production and loss processes are ignorable and denoting the 
ion-electron density by n, the sum of the ion and electron temperatures by T, 
and the flow velocity in the direction of the magnetic field by v, the one dimen- 
sional equations of continuity, momentum and energy for adiabatic flow are: 

n t + v n s + n v s = _ B s (1) 

k T k 

nv t+ nvv s +__n s+ -.nT s = -ng (2) 

n (T t + v T s ) + (1 - y) T (n t + v n s ) = 0. ^ 
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In these equations subscripts denote partial differentiation, B denotes the mag- 
netic field, k is Boltzman's constant, M is the proton mass, g is the magnitude 
of the gravitational acceleration component in the direction of the magnetic field 
and y is the ratio of specific heats. When isothermal flow is considered only 
equations (1) and (2) are applicable with T s = 0. 

The set of hyperbolic equations (1) through (3) are best solved numerically 
by the method of characteristics. Since this is a technique that has been 
thoroughly described in the literature (see for example, Courant and Friedricks, 
1948) the solution scheme will not be described in detail. For the particular 
physical model at hand the numerical solution follows the development employed 
by Grebowsky (1968) and Grebowsky and Smith (1971). 

In order to obtain a complete solution of equations (1) through (3) at all 
times and at all points along a field line, proper boundary conditions must be 
specified. The natural choice for the boundary conditions is a specification of 
n, v and T at all points along the field line at time t = 0 and at all times at the 
base of the field line (i.e., at s = 0). Symmetry of the plasma flow along a trough 
field line about the equatorial plane is assumed so that the condition v = 0 at the 
equator completes the set of required boundary conditions. Due to the assumed 
symmetry of the flow about the magnetic equator the solution will be developed 
explicitly on only one side of the equator. 

The subsonic region on the equator side (i.e., downstream) of the shock is 
coupled directly to the upward directed supersonic flow on the earth side via the 
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shock jump conditions which are derived by requiring the conservation of plasma 
number flux along with the conservation of momentum and energy fluxes through 
the shock discontinuity. Defining the acoustic speed by c 2 = 7 kT/M, the re- 
lationship between the dependent variables on the equator side of the shock (i.e., 
in the subsonic region) denoted by the subscript 1 and those on the earth side of 
the shock, denoted by the subscript 0 , as obtained from the conservation relations 
are: 


n 0 = (y - !) v o + 2 c 0 
n i (7 + 1) V^ 2 


(4) 


n„ 


_i _ 

0 " n i 


(5) 


c l = (— ) [c o - y v o - v')j 


( 6 ) 


The primes denote velocities in the rest frame of the moving shock. That is, 
v’ = v - V s where V s is the shock velocity and v is the observed velocity in the 
frame of reference fixed with respect to the plane of the field line. Using the 
characteristic equations resulting from the hydrodynamic relations ( 1 ) through 
(3) and using relations (4) through (5) to couple the supersonic and subsonic 
regions, the time evaluation of the plasma flow along trough field lines can be 
calculated. 
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Fixed L Coordinate 


Setting the L coordinate at a value of 6, the evaluation of the plasma flow is 
computed from an initial state corresponding to the upward directed supersonic 
steady state flow solution of the hydrodynamic equations from the base of the 
field line to the equator where the shock discontinuity originally forms. The 
nature of similar flow states along field lines characterized by different L 
coordinates can be inferred from these calculations and the previous computa- 
tions of Banks et al. (1971). 

The base level altitude is selected as 3000 kilometers which at trough lati- 
tudes is in the protonosphere (Brintonet al., 1971). Typical values chosen for the 
plasma density and temperature at this altitude where the upward plasma flow 
is assumed to be slightly supersonic are: 

n = 300 ions-electrons/cc 
T = 8000°K 
v = 11.0 km/sec 

The temperature 8000*K corresponds to a sound speed c of 10.49 km/sec for 
adiabatic flow (y = 5/3) and to 8.13 km/sec for isothermal flow (y = 1). These 
values will be used as the assumed time invariant boundary values for all the 
consequent model calculations. 

Model Calculations 

Under the assumption that thermal conduction is sufficient on either side of 
the shock to produce an isothermal state, the flow evolves as shown in Figure 2. 
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The initial upward directed supersonic flow corresponding to steady state condi- 
tions is an upward accelerated flow. As seen in the figure, initially the shock is 
a discontinuity at the equator with vanishing plasma flow and enhanced temper- 
ature and density at the equator. 

As the shock moves away from the equator the supersonic flow on the earth 
side remains unchanged whereas the downstream subsonic flow varies with time. 
The shock velocity decelerates from approximately 12 km/sec when the shock is 
initially at the equator to values less than 0.1 km/sec when the shock reaches 
altitudes near and below 10,000 kilometers. This is in agreement with the study 
of Banks et al., (1971) which found that the shock takes the order of a day or two 
to travel from the equator to the base of the protonosphere. If this model were 
applicable for longer periods of time, the shock would disappear and the flow 
along the entire field line would become subsonic since the decelerating shock 
must become weaker as the flow upstream of the shock decelerates with de- 
creasing altitude towards subsonic values. 

The temperature is enhanced up to a factor of ten fold through the shock 
transition and hence would result in a large scale height for the density variation 
at high altitudes. However, the plasma state departs significantly from diffusive 
equilibrium due to the relatively large magnitude of the subsonic flow velocity 
downstream of the shock. This is particularly evident since the computed vari- 
ation of the subsonic velocity yields an increase of the plasma density with in- 
creasing altitude downstream of the shock. 
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The temperature on the equator side of the shock drops from an initial 
peak of 88,000 e K to values near 24,000°K under isothermal conditions as the 
shock moves earthward. If one arbitrarily assumes that the ion and electron 
temperatures are comparable to one another, the computed temperatures cor- 
respond to individual species energies of the order of a few ev. Energies of this 
order of magnitude have been observed for the thermal ions outside of the 
plasmasphere (Serbu and Maier, 1970). 

If on the other hand the flow was adiabatic along the field line, the flow 
would evolve as shown in Figure 3 from the previously described initial condi- 
tions. It is not likely that the actual flow state is accurately described by the 
adiabatic model since the adiabatic steady state solution is characterized by a 
decrease of the temperature with increasing altitude to values less than 1000'K 
— magnitudes which are much lower than experimentally measured. However, 
the adiabatic model will provide limits for the actual flow and by comparison 
with the isothermal results will contribute to an understanding of the possible 
flow states. 

Again, as for isothermal flow, the shock decelerates as it moves earthward 
and a period of several days is required for the shock to travel from the equator 
to the base of the protonosphere. As the shock moves earthward the density 
everywhere on the equator side increases with time due to the upward plasma 
transport. However, the rate of increase of the ion density at the equator is 
slower than the increase near the shock so that a narrowing depression in the 
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density occurs centered on the equator. Associated with the density depression 
is a peak in the temperature due to the deposition of energy by the equatorially 
converging plasma streams. The adiabatic flow is thus characterized by upward 
directed temperature gradients in the subsonic flow region. 

The previous calculations for a fixed L coordinate indicate the general be- 
havior expected when upward directed supersonic streams form plasma shocks 
on closed magnetic field lines. Both the adiabatic and isothermal models pre- 
dict significant departures from diffusive equilibrium in the expanding subsonic 
region equatorward of the shock. Also both models predict that the shock can 
take several days to travel from the equator to the base of the protonosphere. 

In order to determine whether this behavior is always characteristic of the trough 
flow, the assumption of a time invariant L coordinate must be relaxed in order 
to account for the effects of the more complex convection patterns known to 
exist outside of the plasmasphere. The next section will explore the effects of a 
varying L coordinate. 

Cross L Drift 

The general magnetospheric convection pattern for the thermal plasma out- 
side of the plasmasphere is characterized by cross L drifts. These cross L 
drifts can be explored using the hydrodynamic one dimensional time dependent 
equations to describe the plasma motion along the field line as the field line and 
the plasma frozen to it contract (expand) with a decrease (increase) in the L 
coordinate. 
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The procedures for including the temporal variation of the L coordinate 
within the framework of the one dimensional hydrodynamic equations (1) through 
(3) consists of approximating the continuous change of the L coordinate by a 
discrete variation in which the L coordinate is impulsively changed at small 
time intervals. Assuming that the L coordinate is specified as a function of 
time, the solution scheme for the flow development consists of the following 
steps: 

1. The plasma state is specified at equidistant points along the magnetic 
field line characterized by the coordinate L(t) at time t. With the L 
coordinate fixed at the value L(t) the hydrodynamic and shock relations 
are solved to determine the solution at each computing point at time 

t + At. 

2. The L coordinate is changed to L (t + At) and the dependent variables at 
equidistant points along the L (t + At) line are obtained by a one to one 
mapping of the equidistant computing points along the L(t) line. In this 
transformation the total ion tube content is conserved and the density 
will vary inversely to the incremental change in volume associated with 
each transformed point. For an isothermal state the temperature at a 
point will remain unchanged. However, for adiabatic flow the tempera- 
ture is transformed according to the relation Tn _r+1 = constant. The 
velocity at each point remains unchanged in the transformation except 
at the downstream side of the shock where in order to conserve the total 
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ion content and to preserve the quasistatic nature of the shock process 
the downstream velocity must be computed from the shock relations (4) 
through (6). 

3. Replacing t by t + A t the solution can then be obtained at a later time 
following the above steps. 

Field Line Contraction 

In order to determine what effect a time dependent L coordinate has on the 
field aligned flow, an extreme situation will be treated in which the L variation 
is taken as that characteristic of the inward magnetosphere drift expected at 
midnight for a storm enhanced convective flow state. A simple model often 
used (e.g., Kavanaugh et al., 1968; Chen and Wolf, 1971) to approximate the con- 
vective flow from the magnetosphere tail describes the E x B drift of the ambient 
plasma in the equatorial plane as due to a constant dawn to dusk directed electric 
field. The magnitude of this field is typically of the order of a few tenths of a 
millivolt/meter. Using this model and assuming a dipole field line is character- 
ized initially by the L coordinate L = 11 (i.e., a field line that extends into the 
magnetospheric tail) the variation of L along a convection streamline at midnight 
for E = 0.5 mv/m corresponding to an enhanced flow state is: 


L = 


11 

/2.13 t + 1 


(7) 


This functional variation which ignores the effects of corotation will be employed 
to determine the gross effects of the cross L drift on the development of the high 
speed flow along trough field lines. 
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Again, as previously, the temporal development is computed for both 
adiabatic and isothermal flow states beginning with the steady state solution 
along the L = 11 field line at time t = 0. The isothermal flow development under 
the cross L drift of equation 7 is shown in Figure 4. 

As the plasma convects inward to smaller L coordinates, the flow in the 
supersonic region varies with time as the plasma state continually relaxes 
towards the steady state configurations associated with each contracted field 
line state. Only the velocity variations in the supersonic region are plotted in 
detail in Figure 4 — the corresponding density variations in the supersonic 
region were small compared to the initial profile and were omitted for clarity. 

The temperature on the earth side of the shock was fixed at 8000*K. 

The general features of each of the computed field line states are similar 
to those derived previously with the L coordinate held fixed. The density in the 
subsonic region varies only slightly with altitude along the field lines except in 
the immediate neighborhood of the shock where large gradients in the flow speed 
can develop. The temperature on the equator side of the shock decreases with 
time while the density increases. 

Since the cross L drift model employed approximates the convection inward 
at midnight outside of the plasmasphere, the temporal variations computed at the 
equator correspond to altitude variations. The computed temperature increases with 
altitude attaining values of the order of 90,000 8 K. Also due to the upward trans- 
port of plasma the density decreases with altitude at the equator. This density 


13 



decrease differs significantly from the l/L 4 variation which would exist if the 
cross L drift was not accompanied by a flow of plasma along the field lines. 

The shock, as in previous calculations, decelerates in its motion towards 
the earth. The sensitivity of the shock velocity to variations in the flow resulting 
from the contraction is apparent in the oscillatory variation of the shock velocity. 
Originating at the equator on the L = 11 field line, the shock reaches the base of 
the field line in a period of a few hours at which time it has drifted to an L co- 
ordinate less t han 4. The exact time and location of this transit is of course de- 
pendent upon the functional form of the cross L drift. 

For the same variation of the L coordinate with time an adiabatic flow state 
will evolve as shown in Figure 5. The contraction of the flux tube induces com- 
pression andrarefraction waves which interact with the shock to produce fluctu- 
ations in the magnitude of the shock velocity. After an initial downward acceler- 
ation period near the equator the shock tends to decelerate as it moves into the 
denser plasma near the earth. The shock reaches the base altitude in three 
hours when it has drifted near the L = 4 field line. 

As in the fixed L study the adiabatic flow configurations exhibit more com- 
plex structure than the isothermal development. For 7 = 5/3 the density on the 
equator side of the shock decreases with altitude along the field lines resulting in 
an extensive trough in the density at the equator. The temperature on the other 
hand peaks at the equator as energy is transferred to the equator via the plasma 
motion along the field line and as the plasma is adiabatically heated by the 
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compression of the flux tube. Assuming the field line drifts within a fixed mag- 
netic meridian, the computed temporal development corresponds to an equatorial 
temperature which increases rapidly with decreasing distance from the earth. It 
varies from 15,000°K (of the order of 1 ev) at L = 11 to values near 150,000°K 
(of the order of 10 ev) at L = 4. Hence adiabatic flow produces an extremely hot 
plasma state just outside of the plasmasphere. 

Another interesting feature of the adiabatic development is the production 
of back pressures on the equator side of the shock which drive the subsonic flow 
away from the equator during the initial phase of the shock motion. Eventually 
however the flow throughout the subsonic region is directed upward along the 
field line. Hence in contrast to the isothermal model the subsonic flow can be 
directed in either direction along the field line with the downward flows prevail- 
ing at the lower L coordinates. 

Expanding Flux Tube 

Only a compression of a trough flux tube via convection to smaller L co- 
ordinates was explicitly treated previously. However, the flow model under 
consideration can be extended to include cross L drifts in which the plasma 
convects in the direction of increasing L coordinates. Such drifts occur in the 
magnetosphere since the convective flow from the magnetosphere tail completes 
its cycle by traversing the magnetopause. 

The effects of a drift to higher L coordinates were explored by considering 
the flow evolution when the L coordinate varied according to the reverse of 


15 



equation (7). The initial states for the expansion were chosen to be the solutions 
computed in the previous section at L = 4.8. The results are depicted in Figure 
6. Only the shock velocity and variations at the equator will be considered ex- 
plicitly since the other variations are analogous to those computed previously. 

The initially earthward moving shocks (negative velocities) decelerate until 
the reduction of pressure downstream by the expansion process results in a re- 
versal of the shock motion at L coordinates slightly greater than 5. The shocks 
then accelerate upwards as they drift to higher L coordinates. Thus following a 
typical drift path of a trough field line in the magnetosphere, along which com- 
pression and expansion occur, the cross L drift from the magnetosphere tail 
towards lower L coordinates will be characterized by downward moving shocks; 
whereas, the convection in the direction of increasing L coordinates can be 
characterized by upward moving shocks as the region of supersonic flow expands 
along the field direction towards the steady state polar wind configuration. 

During the expansion under isothermal conditions the equatorial temperature 
initially decreased with increasing L but as the shock moved towards the equator 
the temperature slowly increased. For adiabaticity, on the other hand, the equ- 
atorial temperature always decreased with time since the adiabatic expansion of 
the flux tube cooled the plasma faster than it could be heated by energy transport 
along the field line. In both cases the equatorial density decreased with increas- 
ing distance from the earth. In fact the equatorial density closely follows a l/L 4 
relation indicating that the flux tube expansion affects the density more than the 
transport along the field line. 
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Conclusions 


The evolution of the plasma flow along trough field lines has been computed 
for the case of a fixed L coordinate and for cross L drifts in which the L co- 
ordinate decreases and increases with time. These calculations yield the 
general flow behavior expected when the supersonic trough wind is transformed 
via a shock discontinuity into a subsonic flow. The effect of the shock is to 
produce a region of hot plasma centered about the equator in which the plasma 
can depart significantly from a static diffusive equilibrium state. 

Regions of supersonic flow are always expected along the closed field lines 
with high L coordinates outside of the plasmosphere since the time required for 
the shock to travel from its initial position at the equator to below the base of the 
protonosphere, where the supersonic flow does not exist, is typically much longer 
than the lifetime of a convecting closed field line. The field lines with smaller L 
coordinates near the plasmapause on the other hand may at times be character- 
ized by only subsonic field aligned flows because the shock quickly travels the 
length of these lines of force. Even these subsonic flows however will eventually 
become supersonic as the plasma drifts in the direction of increasing L with an 
attendent upward acceleration which produces supersonic flow and an upward 
moving shock. The expending supersonic region must approach the polar wind 
configuration as the L coordinate increases toward polar values. 

Thus the dynamics of the ambient plasma along trough field lines is seen to 
be very complex with the plasma state along any given field line a function of 
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the past history of its drift motion. Extensive measurements of the ambient 

l 

plasma distribution along the trough field lines are now required to set more 
realistic and flexible boundary conditions and to test these models. 
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Figure 1. Polar wind convects onto closed field lines at night resulting 
in shocks which propagate towards the earth as the plasma drifts around 
the earth to the dayside magnetopause where the polar wind evolves. 
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Figure 2. The temporal development of the plosma flow along an L = 6 magnetic field line is 
plotted for isothermal conditions. The shock is at the equator initially and the supersonic 
flow on the earth side of the shock corresponds to the steady state solution of the hydrody- 
namic equations. 
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Figure 3. The temporal development of the flow along the L = 6 field line is plotted for adiabatic 
flow with the shock beginning its motion at the equator. The supersonic region upstream of the 
shock is the steady state solution. 
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Figure 4. For isothermal flow the field aligned distribution evolves as shown 
when the field line and its frozen in plasma drifts in the direction of increas- 
ing L according to equation 7. The initial distribution along the L = 11 line 
is the steady state supersonic flow. 
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Figure 5. The temporal development of the trough wind under adiabatic conditions is shown for 
the contraction of the flux tube described by equation 7. The initial distribution along the 
L = 11 line is the steady state flow solution with the shock at the equator. 
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Figure 6. The variation of the shock velocity and the equatorial plasma state 
for both adiabatic and isothermal flow is described for a flux tube which ex- 
pands according to the inverse of equation 7. The initial states are the solu- 
tions produced for the contraction of the flux tube to smaller L coordinates. 
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